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COMPLEX LINE BUNDLES OVER SIMPLICIAL COMPLEXES AND 

THEIR APPLICATIONS 

FELIX KNOPPEL AND ULRICH PINKALL 


Abstract. Discrete vector bundles are important in Physics and recently found re¬ 
markable applications in Computer Graphics. This article approaches discrete bundles 
from the viewpoint of Discrete Differential Geometry, including a complete classification 
of discrete vector bundles over finite simplicial complexes. In particular, we obtain a dis¬ 
crete analogue of a theorem of Andre Weil on the classification of hermitian line bundles. 
Moreover, we associate to each discrete hermitian line bundle with curvature a unique 
piecewise-smooth hermitian line bundle of piecewise constant curvature. This is then 
used to define a discrete Dirichlet energy which generalizes the well-known cotangent 
Laplace operator to discrete hermitian line bundles over Euclidean simplicial manifolds 
of arbitrary dimension. 


1. Introduction 

Vector bundles are fundamental objects in Differential Geometry and play an important 
role in Physics [2]. The Physics literature is also the main place where discrete versions 
of vector bundles were studied: First, there is a whole field called Lattice Gauge Theory 
where numerical experiments concerning connections in bundles over discrete spaces (lat¬ 
tices or simplicial complexes) are the main focus. Some of the work that has been done in 
this context is quite close to the kind of problems we are going to investigate here [3, 4, 6]. 

Vector bundles make their most fundamental appearance in Physics in the form of the 
complex line bundle whose sections are the wave functions of a charged particle in a 
magnetic field. Here the bundle comes with a connection whose curvature is given by 
the magnetic field [2]. There are situations where the problem itself suggests a natural 
discretization: The charged particle (electron) may be bound to a certain arrangement of 
atoms. Modelling this situation in such a way that the electron can only occupy a discrete 
set of locations then leads to the “tight binding approximation” [12, 1, 17]. 

Recently vector bundles over discrete spaces also have found striking applications in Ge¬ 
ometry Processing and Gomputer Graphics. We will describe these in detail in Section 2. 

In order to motivate the basic definitions concerning vector bundles over simplicial com¬ 
plexes let us consider a smooth manifold M that comes with smooth triangulation (Fig¬ 
ure 1). 

Let E be a smooth vector bundle over M of rank Then we can define a discrete version 
E of E by restricting E to the vertex set V of the triangulation. Thus E assigns to each 
vertex f G V the .^-dimensional real vector space Ej := Ej. This is the way vector bundles 
over simplicial complexes are defined in general: Such a bundle E assigns to each vertex 
i a .^-dimensional real vector space Ej in such a way that Ej ClEj = 0 for i ^ j- 
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Figure 1. A smooth triangulation of a manifold. 

So far the notion of a discrete vector bundle is completely uninteresting mathematically: 
The obvious definition of an isomorphism between two such bundles E and E just would 
require a vector space isomorphism /*: Ej —)• Ej for each vertex i. Thus, unless we put 
more structure on our bundles, any two vector bundles of the same rank over a simplicial 
complex are isomorphic. 

Suppose now that E comes with a connection V. Then we can use the parallel transport 
along edges ij of the triangulation to define vector space isomorphisms 

• E2 y Ej 

This leads to the standard definition of a connection on a vector bundle over a simplicial 
complex: Such a connection is given by a collection of isomorphisms r]ij : Ej —)• Ej defined 
for each edge ij such that 

Now the classification problem becomes non-trivial because for an isomorphism / between 
two bundles E and E with connection we have to require compatibility with the transport 
maps rfij-. 

fj O flij ~ Vij ° fi- 

Given a connection ry and a closed edge path 7 = • • • ei (compare Section 4) of the 

simplicial complex we can define the monodromy G Aut(Ej) around 7 as 

^7 = o • • • o ^ei • 

In particular the monodromies around triangular faces of the simplicial complex provide 
an analog for the smooth curvature in the discrete setting. In Section 4 we will classify 
vector bundles with connection in terms of their monodromies. 

Let us look at the special case of a rank 2 bundle E that is oriented and comes with a 
Euclidean scalar product. Then the 90°-rotation in each fiber makes it into 1-dimensional 
complex vector space, so we effectively are dealing with a hermitian complex line bundle. 
If ijk is an oriented face of our simplicial complex, the monodromy Pgijk- Ej —)• Ej 
around the triangle ijk is multiplication by a complex number hijk of norm one. Writing 
hijk = with —TT < aijk < vr we see that this monodromy can also be interpreted as a 
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real curvature aijk G (—7r,7r]. It thus becomes apparent that the information provided by 
the connection rj cannot encode any curvature that integrated over a single face is larger 
than ±7r. This can be a serious restriction for applications: We effectively see a cutoff for 
the curvature that can be contained in a single face. 

Remember however our starting point: We asked for structure that can be naturally trans¬ 
ferred from the smooth setting to the discrete one. If we think again about a triangulated 
smooth manifold it is clear that we can associate to each two-dimensional face ijk the 
integral of the curvature 2-form over this face. This is just a discrete 2-form in 
the sense of discrete exterior calculus [5]. Including this discrete curvature 2-form with 
the parallel transport rj brings discrete complex line bundles much closer to their smooth 
counterparts: 

Definition: A hermitian line bundle with curvature over a simplicial complex X is a 
triple (E,r/, n). Here E is complex hermitian line bundle over X, for each edge ij the 
maps r/ij-. Ej —)■ Ej are unitary and the closed real-valued 2-form Q on each face ijk 
satisfies 

Vki o rijk o Vij = idp,- 

In Section 7 we will prove for hermitian line bundles with curvature the discrete analog 
of a well-known theorem by Andre Weil on the classification of hermitian line bundles. 

In Section 8 we will define for hermitian line bundles with curvature a degree (which can 
be an arbitrary integer) and we will prove a discrete version of the Poincare-Hopf index 
theorem concerning the number of zeros of a section (counted with sign and multiplicity). 

Finally we will construct in Section 10 for each hermitian line bundle with curvature a 
piecewise-smooth bundle with a curvature 2-form that is constant on each face. Sections 
of the discrete bundle can be canonically extended to sections of the piecewise-smooth 
bundle. This construction will provide us with finite elements for bundle sections and 
thus will allow us to compute the Dirichlet energy on the space of sections. 

2. Applications of Vector Bundles in Geometry Processing 

Several important tasks in Geometry Processing (see the examples below) lead to the 
problem of coming up with an optimal normalized section f of some Euclidean vector 
bundle E over a compact manifold with boundary M. Here “normalized section” means 
that (f is defined away from a certain singular set and where defined it satishes \(f\ = 1. 

In all the mentioned situations E comes with a natural metric connection V and it turns 
out that the following method for finding f yields surprisingly good results: 

Among all sections V' o/ E find one which minimizes /jvi under the constraint 

Jm IV'P = 1- Then away from the zero set of use 4> = V'/IV’I- 

The term ’’optimal” suggests that there is a variational functional which is minimized by 
(p and this is in fact the case. Moreover, in each of the applications there are heuristic 
arguments indicating that (p is indeed a good choice for the problem at hand. For the 
details we refer to the original papers. Here we are only concerned with the Discrete 
Differential Geometry involved in the discretization of the above variational problem. 

2.1. Direction Fields on Surfaces. Here M is a surface with a Riemannian metric, E = 
TM is the tangent bundle and V is the Levi-Civita connection. Figure 2 shows the resulting 
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unit vector field If we consider TM as a complex line bundle, normalized sections of the 



Figure 2. An optimal direction field on a surface. 

tensor square L = TM (8) TM describe unoriented direction fields on M. Similarly, “higher 
order direction fields” like cross fields are related to higher tensor powers of TM. Higher 
order direction fields also have important applications in Computer Graphics. 

2.2. Stripe Patterns on Snrfaces. A stripe pattern on a surface M is a map which 
away from a certain singular set assigns to each point p G M an element (j){p) G S = {z G 
C|| 2 ;| = 1}. Such a map cj) can be used to color M in a periodic fashion according to a 
color map that assigns a color to each point on the unit circle S. Suppose we are given a 
1-form w on M that specihes a desired direction and spacing of the stripes, which means 
that ideally we would wish for something like cj) = e*" with da = uo. Then the algorithm in 
[9] says that we should use a cj) that comes from taking E as the trivial bundle E = M x C 
and Vip = dll’ — iujip. Sometimes the original data come from an unoriented direction field 
and (in order to obtain the 1-form ca) we first have to move from M to a double branched 
cover M of M. This is for example the case in Figure 3. 



Figure 3. An optimal stripe pattern aligned to an unoriented direction field. 
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2.3. Decomposing Velocity Fields into Fields Generated by Vortex Filaments. 

The velocity fields that arise in fluid simulations quite often can be understood as a 
superposition of interacting vortex rings. It is therefore desirable to have an algorithm 
that reconstructs the underlying vortex filaments from a given velocity held. Let the 
velocity held v on a domain M C be given as a 1-form oj = (v, •). Then the algorithm 
proposed in [20] uses the function (/>: M —)• C that results from taking the trivial bundle 
E = M X C endowed with the connection Vi/^ = dip — iujip. Note that so far this is just 
a three-dimensional version of the situation in Section 2.2. This time however we even 
forget (p in the end and only retain the zero set of ip as the hlament conhguration we are 
looking for. 



Figure 4. A knotted vortex hlament dehned as the zero set of a complex 
valued function ip. It is shown as the intersection of the zero set of Help 
with the zero set of Ivciip. 


2.4. Close-To-Conformal Deformations of Volumes. Here the data are a domain 
M C and a function n: M —)■ M. The task is to hnd a map /: M —>■ which is 

approximately conformal with conformal factor e“, i.e. for all tangent vectors X G TM 
we want 

\df{X)\^e-\X\. 

The only exact solutions of this equations are the Mobius transformations. For these we 
hnd 

df{X) = e^PpXiP 

for some map "0: M —)• HI with |^/;| = 1 which in addition satishes 

dip{X) = — 2 (gradrt x X) ip. 

Note that here we have identihed with the space of purely imaginary quaternions. Let 
us dehne a connection V on the trivial rank 4 vector bundle M x H by 

VxV’ := dip{X) -|- |(gradrt x X)ip. 

Then we can apply the usual method and hnd a section (/>: M —)• BI with \(p\ = 1. In 
general there will not be any /: M —)• that satishes 

(2.1) df{X) = 

exactly but we can always look for an / that satishes (2.1) in the least squares sense. See 
Figure 5 for an example. 
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Figure 5. Close-to-conformal deformation of a sphere based on a desired 
conformal factor specified as the potential of a collection of point charges. 


3. Discrete Vector Bundles with Connection 

An (abstract) simplicial complex is a collection X of finite non-empty sets such that if a 
is an element of X so is every non-empty subset of a ([15]). 

An element of a simplicial complex X is called a simplex and each non-empty subset of 
a simplex a is called a face of a. The elements of a simplex are called vertices and the 
union of all vertices V = Uo-ex t is called the vertex set of X. The dimension of a simplex 
is defined to be one less than the number of its vertices: dimer := |(t| — 1. A simplex of 
dimension k is also called a k-simplex. The dimension of a simplicial complex is defined 
as the maximal dimension of its simplices. 

To avoid technical difficulties, we will restrict our attention to finite simplicial complexes 
only. The main concepts are already present in the finite case. Though, the definitions 
carry over verbatim to infinite simplicial complexes. 

Definition 1. Let ¥ be a field and let X be a simplicial complex with vertex set V. A 
discrete ¥-vector bundle of rank G N over X is a map tt: E —)• V such that for each 
vertex i G V the fiber over i 

E. :=vr-H{i}) 

has the structure of a A-dimensional ¥-vector space. 

Most of the time, we slightly abuse notation and denote a discrete vector bundle over a 
simplicial complex schematically by E —)• X. 

The usual vector space constructions carry fiberwise over to discrete vector bundles. So we 
can speak of tensor products or dual bundles. Moreover, the fibers can be equipped with 
additional structures. In particular, a real vector bundle whose fibers are Euclidean vector 
spaces is called a discrete Euclidean vector bundle. Similarly, a complex vector bundle 
whose fibers are hermitian vector spaces is called a discrete hermitian vector bundle. 

So far discrete vector bundles are completely uninteresting from the mathematical view¬ 
point - the obvious definition of an isomorphism / between two discrete vector bundles 
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E and E just requires isomorphisms between the fibers /i: Ej —)■ E*. Thus any two bun¬ 
dles of rank ^ are isomorphic. This changes if we connect the fibers along the edges by 
isomorphisms. 

Let a = {io,... ,ik} be a A:-simplex. We define two orderings of its vertices to be equiv¬ 
alent if they differ by an even permutation. Such an equivalence class is then called an 
orientation of a and a simplex together with an orientation is called an oriented simplex. 
We will denote the oriented fc-simplex just by the word io---ik- Eurther, an oriented 
1-simplex is called an edge. 

Definition 2. Let E —)■ X 6e a discrete vector bundle over a simplicial eomplex. A diserete 
conneetion on Li is a map rj which assigns to each edge ij an isomorphism rjij ; Ej —)■ Lj 
of vector spaces sueh that 

hji = 

Remark 1: Here and in the following a morphism ofveetor spaces is a linear map that also 
preserves all additional structures - if any present. E.g., if we are dealing with hermitian 
veetor spaces, then a morphism is a eomplex-linear map that preserves the hermitian 
metric, i.e. it is a complex linear isometric immersion. 

Definition 3. A morphism of discrete vector bundles with eonnection is a map / : E —)■ E 
between discrete vector bundles E —)■ X and E —)■ X with eonnections rj and 6 (resp.) such 
that 

i) for each vertex i we have that /(Ej) C Ej and the map fi = /l^v Ej —)■ Ej is o 
morphism of vector spaces, 
a) for each edge ij the following diagram commutes: 



i.e. 9ij o fi = fj o riij. 

An isomorphism is a morphism which has an inverse map, which is also a morphism. 
Two discrete vector bundles with connection are called isomorphic, if there exists an 
isomorphism between them. 

Again let V denote the vertex set of X. A discrete vector bundle E —)• X with connection 
T] is called trivial, if it is isomorphic to the product bundle 

;= V X F-^ 

over X equipped with the connection which assigns to each edge the identity id^j?. 

It is a natural question to ask how many non-isomorphic discrete vector bundles with 
connection exist on a given simplicial complex X. 

Remark 2: For /c G N, the k-skeleton X^ of a simplicial complex X is the subcomplex that 
consists of all simplices a £ X of dimension < k, 

X^ := jcj G X I dim O' < fej. 

The classification of vector bundles over X only involves its 1-skeleton X^ and could be 
equally done just for discrete veetor bundles over 1-dimensional simplicial complexes, i.e. 
graphs. Later on, when we consider discrete hermitian line bundles with curvature in 
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Section 1, the 2- and H-skeleton come into play and finally, in Section 11, we will use the 
whole simplicial complex. 


4. Monodromy - A Discrete Analogue of Kobayashi’s Theorem 

Let X be a simplicial complex. Each edge ij of X has a start vertex s{ij) := i and a 
target vertex t{ij) := j. A edge path 7 is a sequence of successive edges (ei,..., efij, i.e. 
s(efc+i) = t{ek) for all A: = 1, 1, and will be denoted by the word: 

7 = e£ - • - ei. 

If i = s(ei), we say that 7 starts at i, and if j = t{ei) that 7 ends at j. The complex X 
is called connected, if any two of its vertices can be joined by an edge path. From now on 
we will only consider connected simplicial complexes. 

Now, let E —7- X be a discrete vector bundle with connection 7. To each edge path 
7 = • • • ei from i to j, we define the parallel transport Ey: E* —)■ Ej along 7 by 

P-, := r]ee ° ° hei ■ 

To each edge path 7 = • • • ei we can assign an inverse path 7“^. If 7 = Cm • • • e^+i starts 

where 7 ends, we can build the concatenation 77: With the notation ij~^ := ji, we have 

7"^ := • • • ej^, 77 := • • • e^+ie^ • • • ei. 

Whenever 77 is defined, 

(4.1) = P ^ o P ^, P^-i = Pfi ^■ 

The elements of the fundamental group are identified with equivalence classes of edge 
loops, i.e. edge paths starting and ending a given base vertex i of X, where two such loops 
are identified if they differ by a sequence of elementary moves [?]: 

ec ■ ■ ek+ie~^e • • • ei <—)• • • • e^+iefc • • • ei. 

Now, by Equation (4.1), we see that the parallel transport descends to a representation 
of the fundamental group 7ri(X^,i). We encapsulate this in the following 

Proposition 1. Let E —)■ X 6e a discrete vector bundle with connection over a connected 
simplicial complex. The parallel transport descends to a representation of the fundamental 
group 7ri(X^,i).- 

911: 7ri(X^,i) —)■ Aut(Ej), [7] 1—)■ Py. 

The representation 911 will be called the monodromy of the discrete vector bundle E. 

If we change the base vertex this leads to an isomorphic representation - an isomorphism 
is just given by the parallel transport Py along an edge path joining the two base vertices. 
Moreover, if /: E —)■ E is an isomorphism of discrete vector bundles with connection over 
the simplicial complex X, for any edge path 7 = • • • ei from a vertex i to a vertex j the 

following equality holds: 

Py = fjOP,of-\ 

Here P and P denote the parallel transports of E and E. Thus we obtain: 

Proposition 2. Isomorphic discrete vector bundles with connection have isomorphic mon- 
odromies. 
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In fact, the monodromy completely determines a discrete vector bundle with connection 
up to isomorphism, which provides a complete classihcation of discrete vector bundles 
with connection: Let X be a connected simplicial complex. Let E —)■ X be a discrete 
F-vector bundle of rank with connection and let 9JI: 7 ri(X^,i) —)■ Aut(Ej) denote its 
monodromy. Any choice of a basis of the fiber Ej determines a group homomorphism p G 
Hom( 7 ri(X^, i), GL(.^, F)). Any other choice of basis determines a group homomorphism 
p which is related to p by conjugation, i.e. there is 5 G GL(.^, F) such that 

^([ 7 ]) = S • p([7]) • S~^ for all [7] G 7ri(X\i). 

Hence the monodromy XU determines a well-defined conjugacy class of group homomor- 
phisms from 7 ri(X^,i) to GL(.^, F), which we will denote by [XU]. The group GL(.^,F) will 
be referred to as the structure group of E. 

Let XJf(X) denote the set of isomorphism elasses of ¥-vector bundles of rank ^ with 
conneetion overX and let Hom( 7 ri(X^, i), GL(.^, F))/..., denote the set of eonjugacy classes 
of group homomorphisms from the fundamental group 7ri(X^,i) into the strueture group 
GL(j^,F). 

Theorem 1. F: XJip(X) —)■ Hom(7ri(X^, i), GL(.^, F))/...., [E] 1 —)■ [XR] is bijective. 


Proof. By Proposition 2, F is well-defined. Eirst we show injectivity. Gonsider two discrete 
vector bundles E, E over X with connections 7 , fj and let XR, XR denote their monodromies. 
Suppose that [XU] = [fOT]. If we choose bases {Xf, ..., X^} of Ej and {Xl, ..., Xf} of Ej, 
then XU and iOT are represented by group homomorphisms p, p £ Hom( 7 ri(X^, i), GL(.ft, F)) 
which are related by conjugation. Without loss of generality, we can assume that p = p- 
Now, let T be a spanning tree of X with root i. Then for each vertex j of X there 
is an edge path yjj from the root i to the vertex j entirely contained in T. Since T 
contains no loops, the path 7jj is essentially unique, i.e. any two such paths differ by 
a sequence of elementary moves. Thus we can use the parallel transport to obtain bases 
{Xj, ..., X^} C Fj and {X^, • • •, X^} C Fj at every vertex j of X. With respect to these 
bases the connections 7 and rj are represented by elements of GL(.^, F). By construction, 
for each edge e in T the connection is represented by the identity matrix. Moreover, to 
each edge e = jk not contained in T there corresponds a unique [ 7 e] G 7 ri(X^, f). With the 
notation above, it is given by je = if 1^1 ip- particular, on the edge e both connections 
are represented by the same matrix p([ 7 e]) = p{[le\)- Thus, if we define / : E —)■ E by 
f{X'l) := A™ for m = 1,...,.^, we obtain an isomorphism, i.e. E = E. Hence F is 
injective. 

To see that F is surjective we use T to equip the product bundle E := V x F-*^ with a 
particular connection 7. Let p G Hom( 7 ri(X^,z), GL(.^,F)). If e lies in T we set 7 e = id 
else we set pe '■= p{[le\)- By construction, E([E]) = [p]. Thus F is surjective. □ 

Remark 3: Note that Theorem 1 can be regarded as a discrete analogue of a theorem of S. 
Kobayashi f'[10, 14]J, whieh states that the equivalence elasses of connections on principal 
G-bundles over a manifold M are in one-to-one correspondence with the conjugacy classes 
of eontinuous homomorphisms from the path group 4>(M) to the structure group G. In fact, 
the fundamental group of the 1-skeleton is a diserete analogue o/<I>(M). 
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5. Discrete Line Bundles - The Abelian Case 

In this section we want to focus on discrete line bundles, i.e. discrete vector bundle of rank 
1. Here the monodromy descends to a group homomorphism from the closed 1-chains to 
the multiplicative group F* := F \ {0} of the underlying field. This leads to a description 
by discrete differential forms (Section 6). 

Let L —7- X be discrete F-line bundle over a connected simplicial complex. In this case the 
structure group is just F*, which is abelian. Thus we obtain 

Hom(7ri(X\i),F*= Hom(^7ri(X^, i),F*). 

Hom(7ri(X^, i), F*) carries a natural group structure. Moreover, the isomorphism classes 
of discrete line bundles over X form an abelian group. The group structure is given by 
the tensor product: For [L], [L] G QJjp(X), we have 

[L][L] = [L®L], [L]-^ = [L*]. 

The identity element is given by the trivial bundle. In the following we will denote the 
group of isomorphism classes of ¥-line bundles over X by 

The map F: —)■ Hom(^7ri(X^, i),F*), [L] i—)■ [911] is a group homomorphism. By Theo¬ 

rem 1, T is then an isomorphism. 

Now, since F* is abelian, each homomorphism p G Hom(^7ri(X^, i),F*) factors through the 
abelianization 

TTi{X^,i)ab = 7ri(X\i)/[7ri(X\i),7ri(X^i)], 

i.e. for each p G Hom(7ri(X^,i),F*) there is a unique pab £ Hom(7ri(X^,ijab,F*) such 
that 

P ~ Pab ° '^ab- 

Here iTab'- 7ri(X^,i)afe denotes the canonical projection. This yields an isomor¬ 

phism between Hom(7ri(X^, i),F*) and Hom(7ri(X^, i)ab,F*). In particular, 

Tx - Hom(7ri(X\i)ab,F*). 

As we will see below, the abelianization 7ri(X^,i)afe is naturally isomorphic to the group 
of closed 1-chains. 

The group of k-chains Cfc(X,Z) is defined as the free abelian group which is generated by 
the /c-simplices of X. More precisely, let denote the set of oriented k-simplices of X. 
Clearly, for k > 0, each /c-simplex has two orientations. Interchanging these orientations 
yields a fixed-point-free involution pk : Xff —)■ X^. The group of A:-chains is then explicitly 
given as follows: 

Cfc(X,Z) := {c: X^^ ^ Z | co pk = -c). 

Since simplices of dimension zero have only one orientation, Xq’’ = Xq. Thus, 

Co(X,Z) := {c: X^’' ^Z}. 

It is common to identify an oriented A:-simplex a with its elementary k-chain, i.e. the 
chain which is 1 for a, —1 for the oppositely oriented simplex and zero else. With this 
identification a A:-chain c can be written as a formal sum of oriented /c-simplices with 
integer coefficients: 

m 

c = Y^ niai, ni gZ, ai G X°{. 

i=l 
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The boundary operator d^: Cfc(X,Z) —)• Cfc_i(X,Z) is then the homomorphism which is 
uniquely determined by 

k 

dk io ' ' ' ik — ^ ' ( 1 )'^ ^0 ■ ■ ■ • 

j=0 

It well-known that dk o dk+i = 0. Thus we get a chain complex 

0 A Co(x,z) A Ci(x,z) A ••• A Cfc(x,z) A±^ •••. 

The simplicial Homology groups Hfc(X,Z) may be regarded as a measure for the deviation 
of exactness: 

Hfc(X,Z) := kerdk/imdk+i- 

The elements of ker are called k-cycles, those of im9fc_|_i are called k-boundaries. 

It is a well-known fact that the abelianization of the first fundamental group is the first 
homology group (see e.g. [7]). Now, since the first homology of the 1-skeleton consists 
exactly of all closed chains of X, we obtain 

7 ri(X\i)afe = kercli. 

The isomorphism is induced by the map 7 ri(X^,f) —)■ kerSi given by [ 7 ] 1 —)■ where 

7 = • • • ei- We summarize the above discussion in the following theorem. 

Theorem 2. The group of isomorphism classes of line bundles is naturally isomorphic 
to the group Hom(ker 5i,F*).- 

= Hom(ker 9i,F*). 

The isomorphism of Theorem 2 can be made explicit using discrete F*-valued 1-forms 
associated to the connection of a discrete line bundle. 

6. Discrete Connection Forms 
Throughout this section X denotes a connected simplicial complex. 

Definition 4. Let 6 be an abelian group. The group of (5-valued discrete k-forms is 
defined as follows: 

D^(X, 6 ) := |cj: Cfc(X) —)■ 6 | a; group homomorphism^. 

The discrete exterior derivative dk is then defined to be the adjoint of dk+i, i.e. 

dk : D^(X, 0) —^ ^^■'■^(X, 0), dkio := io o dk+i- 

By construction, we immediately get that dk+i o = 0. The corresponding cochain 
complex is called the discrete de Rahm complex with coefficients in 0: 

0 ^ D°(X, 0) A D^(X, 0) A • • • A 4 6 ) A • • • . 

In analogy to the construction of the homology groups, the k-th de Rahm Cohomology 
group H^(X, 0) with coefficients in 0 is defined as the quotient group 

H^(X, 0) := ker dfc/imdfc_i. 

The discrete /c-forms in kerd^ are called closed, those in irnd^-i are called exact. 

Now, let 0L denote the space of connections on the discrete F-line bundle L —)■ X: 

■■= {7 I 7 connection on L| . 
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Any two connections r],9 £ € 1 , differ by a unique discrete 1-form a; G r2^(X,F*); 

9 = Lorf. 

Hence the group F*) acts simply transitively on the space of connections (Hl- Ih 

particular, each choice of a base connection /3 G £l establishes an identification 

€-1,3 rj = ujjd i —> oj G n^(X,F*). 

Remark 4: Note that each discrete vector bundle admits a trivial connection. To see this 
choose for each vertex a basis of the corresponding fiber. The corresponding coordinates 
establish an identification with the product bundle. Then there is a unique connection that 
makes the diagrams over all edges commute. 

Definition 5. Let rj G £l- A connection form representing the connection r] is a 1-form 
OJ G n^(X,F*) such that rj = ojfi for some trivial base connection (3. 

Clearly, there are many connection forms representing a connection. We want to see how 
two such forms are related. 

More generally, two connections ij and 9 in € 1 , lead to isomorphic discrete line bundles if 
and only if for each fiber there is a vector space isomorphism /*: Lj —Lj, such that for 
each edge ij: 

^ij ^ fi ~ fj ° Vij- 

Since pe and 9e are linear, this boils down to discrete F*-valued functions and the relation 
characterizing an isomorphism becomes 

0ij = {djdi )'nij — idg)ijVijj 

i.e. Ij and 9 differ by an exact discrete F*-valued 1-form. In particular, the difference of 
two connection forms representing the same connection rj is exact. 

Thus we obtain a well-defined map sending a discrete line bundle L with connection to 
the corresponding equivalence class of connection forms 

M G n\X,¥^)/dkf{X,¥^). 

Theorem 3. The map F: -3 H^(X, F*)/(in^(X,F*), [L] 1 —)■ [oj\, where u is a connec¬ 

tion form of L, is an isomorphism of groups. 

Proof. Clearly, F is well-defined. Let L and L be two discrete complex line bundle with 
connections rj and 9, respectively. If /? G CIl and f3 G are trivial, so is /I (8) /3 G 
Hence, with rj = ojfi and fj = ujj3, we get 

F([L® L]) = M = M[cD] = F([L])F([L]). 

By the preceding discussion, F is injective. Surjectivity is also easily checked. □ 

Next we want to prove that the map given by 

L!^(X,F*)/dfl°(X,F*) 9 [cj] ^ wikerai ^ Hom(kerai,F*) 

is a group isomorphism. Clearly, it is a well-defined group homomorphism. We show its 
bijectivity in two steps. First, the surjectivity is provided by the following 

Lemma 1. Let X be a simplicial complex and <5 be an abelian group. Then the restriction 
map <I>: 11^(X, 6) —Hom(ker 9^, 0), oo 1 —)• is surjective. 
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Proof. If we choose an orientation for each simplex in X, then is given by an integer 
matrix. Now, there is a unimodular matrix U such that d^U = {0\H) has Hermite 
normal form. Write U = {A\B), where d^A = 0 and d^B = H and let a, denote the 
columns of A, i.e. A = (ai,...,a£). Clearly, a* G ker5fc. Moreover, if c G kerdk, then 
0 = dkC = {0\H)U~^c. Hence U~^c = {q, 0)"'', q G Z^, and thus c = Aq. Therefore {oj | i = 
is a basis of ker d^- Now, let /x G Hom(ker<9^,Z). A homomorphism is completely 
determined by its values on a basis. We define u = (;u(ai),..., fr{ae),0 ..., 0)U~^. Then 
uj G H^(X,Z) and uoA = (^(oi),... ,/x(a^)). Hence <h(w) = /U and $ is surjective for forms 
with coefficients in Z. Now, let 0 be an arbitrary abelian group. And fx G Hom(ker 5^, 0). 
Now, if oi,.., ai is an arbitrary basis of kercifc, then there are forms wi, ... ,uii G H^(X,Z) 
such that 0 Ji{aj) = <5^. Since Z acts on 0, we can multiply Ui with elements G 0 to 
obtain forms with coefficients in 0. Now, set oj = Y^i=iOJi ■ Then lo G H^(X, 0) 

and uj{ai) = fx{ai) ioi i = 1,... ,i. Thus <h(w) = fx. Hence <I> is surjective for forms with 
coefficients in arbitrary abelian groups. □ 

The injectivity is actually easy to see; If = 0, we define an F*-valued function / 

by integration along paths: Fix some vertex i. Then 



where 7 is some path joining i to j. Since ~ 0, the value /(j) does not depend 

on the choice of the path 7. Moreover, df = ux. Together with Lemma 1, this yields the 
following theorem. 

Theorem 4. The map F: H^(X,F*)/dH°(X,F*) — Hom(ker(9i,F*), [w] 1 —)■ is an 

isomorphism of groups. 

Now, let us relate this to Theorem 2. Let L —)■ X be a line bundle with connection ij and 
w be a connection form representing rj, i.e. r] = ojjd for some trivial base connection /3. 
Let [7] G 7ri(X^,i), where 7 = By linearity and since trivial connections have 

vanishing monodromy, we obtain 

9^([7]) = ??e<. o • • • o hei = • • • Wei • /3e^ O • • • O = ^(^^^([y])) • MIl, . 

Hence, by the uniqueness of [9II]ab, we obtain the following theorem that brings everything 
nicely together. 

Theorem 5. Let L —)■ X 6e a line bundle with connection p. Let 9JI denote its monodromy 
and let oj he some connection form representing rj. Then, with the identifications above, 

[mu = [oj]. 

7 . Curvature - A Discrete Analogue of Weil’s Theorem 

In this section we describe complex and hermitian line bundles by their curvature. For 
the first time we use more than the 1-skeleton. 

Let X be a connected simplicial complex and 0 an abelian group. Since = 0, the 
exterior derivative descends to a well-defined map on H^(X, 0)/dH^“^(X, 0), which again 
will be denoted by d. Explicitly, 

d: L!^(X,0)/dH'=-i(X,0) ^ L!^+^(X,0), [w] ^ dux. 
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Definition 6 . The ¥ ^-curvature of a discrete ¥-line bundle L —)■ X is the discrete 2-form 
ri G rj^(X, F*) given by 

Q = d[uj], 

where [w] G r2^(X,F*)/(iO°(X, F*) represents the isomorphism class [L]. 

Remark 5: Note that Q just encodes the parallel transport along the boundary of the 
oriented 2-simplices of X - the “local monodromy”. 

From the definition it is obvious that the F*-curvature is invariant under isomorphisms. 
Thus, given a prescribed 2-form Q G r2^(X,F*), it is a natural question to ask how many 
non-isomorphic line bundles have curvature 

Actually, this question is answered easily: If d[u}] =17 = d[Cj], then the difference of u and 
oj is closed. Factoring out the exact 1-forms we see that the space of non-isomorphic line 
bundles with curvature Q can be parameterized by the first cohomology group H^(X,F*). 
Furthermore, the existence of a line bundle with curvature 17 G 17^(X,F*) is equivalent to 
the exactness of 17. 

But when is a /c-form 17 exact? Certainly it must be closed. Even more, it must vanish 
on every closed fc-chain: If 17 = im d and S' is a closed fc-chain, then 

17(5) = dcoiS) = uj(dS) = 0. 

For /c = 1, as we have seen, this criterion is sufficient for exactness. For k > 1 this is not 
true with coefficients in arbitrary groups. 

Example: Consider a triangulation X of the real projective plane MP^. The zero-chain is 
the only closed 2-chain and hence each 'Z 2 -valued 2-form vanishes on every closed 2-chain. 
But H^(X,Z2) = Z2 and hence there exists a non-exact 2-form. 

In the following we will see that this cannot happen for fields of characteristic zero or, 
more generally, for groups that arise as the image of such fields. 

Clearly, there is a natural pairing of Z-modules between 17^(X, 0) and Cfe(X,Z); 

(.,.): n^{X, (3) X Cfc(X, Z) ^ 6, (w, c) ^ a;(c). 

This pairing is degenerate if and only if all elements of 0 have bounded order. In partic¬ 
ular, if 0 is a field F of characteristic zero, (.,.) yields a group homomorphism 

Fk-. Cfc(X,Z) ^ HomF(17*^(X,F),F) = (17^(X,F))T 

A basis of Cfc(X,Z) is mapped under F^ to a basis of (17^(X,F))* and hence Cfe(X,Z) 
appears as an n^-dimensional lattice in (17^(X,F))*. 

Let df denote the adjoint of the discrete exterior derivative dk with respect to the natural 
pairing between 17^(X,F) and (17^(X,F))*. Clearly, 

dl o Fk = Fk o dk+i. 

Now, since the simplicial complex is finite, we can choose bases of Cfc(X,Z) for all k. This 
in turn yields bases of (I7^(X, F))* and hence, by duality, bases of f7^(X,F). With respect 
to these bases we have 

(7.1) Cfc(X,Z) = Z’"'' C F’"'' = (I7^(X,F))* = I7^(X,F), 

where Uk denotes the number of /c-simplices. Moreover, the pairing is represented by the 
standard product. The operator d’l_^ = dk is then just an integer matrix and 

dk—\. 
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Figure 6. With the identifications 7.1, the space of /c-forms becomes a 
direct sum of the image of dk-i and the kernel of its adjoint the latter 
of which contains the closed /c-chains as a lattice. 

We have imdk-i -L kerd^_j^. Moreover, by the rank-nullity theorem, 

Uk = dimimd^_;^ -|-dimker = dimimdfc-i -|-dimker 

Hence, under the identifications above, we have that F”* = imdk-i ©ker(i|,_^ (see Fig¬ 
ure 6). Moreover, ker contains a basis of ker(i^_^. From this we conclude immediately 
the following lemma. 

Lemma 2. Let oj G Ll^{X,¥), where ¥ is a field of characteristic zero. Then 
uj G imdfc-i (w,c) = 0 for all c G kercf^. 

Remark 6: Note, that for boundary cycles the condition is nothing hut the closedness of 
the form u. Thus Lemma 2 states that a closed form oj G Ll^{X,¥) is exact if and only if 
the integral over all homology classes [c] G Hfc(X, Z) vanishes. 

Let 6 be an abelian group. The sequence below will be referred to as the k-th fundamental 
sequence of forms with coefficients in ©: 

©) n'^iX, ©) % Hom(ker5fc, ©) ^ 0, 
where denotes the restriction to the kernel of dk, i.e. d>fc(w) ;= 

Combining Lemma 1 and Lemma 2 we obtain that the fundamental sequence with coef¬ 
ficients in a field F of characteristic zero is exact for all /c > 1. This serves as an anchor 
point. The exactness propagates under surjective group homomorphisms. 

f 

Lemma 3. Let 21 —> —)■ 0 be a an exact sequence of abelian groups. Then, if the k-th 
fundamental sequence of forms is exact with coefficients in 21, so it is with coefficients in 

25. 

Proof. By Lemma 1 the restriction map is surjective for every abelian group. It is 
left to check that ker<hfc = irnd^-i with coefficients in 25. Let H G 25) such that 

= 0. Since /: 21 —)■ 25 is surjective, there is a form S G Ll^{X, 21) such that = /oS. 
Since 0 = <hfc(rf) = / o <kfc(S), we obtain that ^fc(S) takes its values in ker/. Since is 
surjective for arbitrary groups, there is 0 G Ll^{X,'ker f ) such that <hfc(S) = $^(0). Hence 
— 0) = 0. Thus there is a form ^ G 21) such that dk-if = S — 0. Now, let 
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00 := fo^e 0^-1 (X, 53). Then 

dk-iu: = dk-if o ^ = / o dk-iC = fo{E-G) = foE = n. 

Hence ker<hfc = imdk-i and the sequence (with coefficients in 53) is exact. □ 

Remark 7: The map / ; C —)• C, z i—)• exp(27^^2:) provides a surjective group homomor¬ 
phism from C onto C*, and similarly from M onto S. Hence the k-th fundamental sequence 
of forms is exact for coefficients in C* and in the unit circle S. 

Remark 8: The k-th fundamental sequence with coefficients in an abelian group 0 is exact 
if and only if Of {%,&)/= Honi(ker 5^., 0). The isomorphism is induced by 
the restriction map 

The following corollary is a consequence of Remark 7. It nicely displays the fibration of 
the complex line bundles by their C*-curvature. 

Corollary 1. For 0 = S, C* the following sequence is exact: 

1 ^ h1(X,0) --a lli(X,0)/dll°(X,0) 4 fl2(X,0) ^ Hom(kera2,0) ^ 1. 

Definition 7. Let G H^(X,S). A real-valued form Ll G H^(X,M) is called compatible 
with Ll* if LT = exp(zH). A discrete hermitian line bundle with curvature is a discrete 
hermitian line bundle L with connection equipped with a closed 2-form compatible with the 
S-curvature o/L. 

For real-valued forms it is common to denote the natural pairing with the /c-chains by an 
integral sign, i.e. for iv G n^(X, M) and c G Cfc(X, Z) we write 

/ cj := (w, c) = a;(c). 

J C 

Theorem 6. Let L be a discrete hermitian line bundle with curvature Q. Then Q is 
integral, i.e. 

[ n G 27rZ, for all c G kercl2. 


Proof. By definition the curvature form H satisfies exp(^^l) = duj for some connection 
form (jj G n^(X, S). Thus, if c G ker52, 


exp(z/ fl) = (exp(zn),c) = {duj,c) 


{uj, dc) = 1. 


This proves the claim. 


□ 


Conversely, Corollary 1 yields a discrete version of a theorem of Andre Weil ([19, 11]), 
which states that any closed smooth integral 2-form on a manifold M can be realized as 
the curvature of a hermitian line bundle. This plays a prominent role in the process of 
prequantization [18]. 

Theorem 7. If 11^(X,M) is integral, then there exists a hermitian line bundle with 
curvature Ll. 


Proof. Consider LI* := exp{iLl). Since LI is integral, {Ll*,c) = 1 for all c G ker92- By 
Corollary 1, there exists r G n^(X,S) such that dr = 11*. This in turn defines a hermitian 
line bundle with curvature H. □ 
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Remark 9: Moreover, Corollary 1 shows that the connections of two such bundles differ 
by an element o/H^(X, S). Thus the space of discrete hermitian line bundles with fixed 
curvature can be parameterized by H^(X,S). 

8. The Index Formula for Hermitian Line Bundles 

Before we define the degree of a discrete hermitian line bundle with curvature or the 
index form of a section, let us first recall the situation in the smooth setting. Therefore, 
let L —7- M be a smooth hermitian line bundle with connection. Since the curvature 
tensor of V is a 2-form taking values in the skew-symmetric endomorphisms of L, it 
is completely described by a closed real-valued 2-form H G 

= -zH. 

The following theorem shows an interesting relation between the index sum of a section 
if G r(L), the curvature 2-form H, and the rotation form ofif. This form is defined as 
follows; 

{ip,if) 

Theorem 8. Let L —)■ M 6e a smooth hermitian line bundle with connection and H its 
curvature 2-form. Let if G r(L) be a section with a discrete zero set Z. Then, if C is a 
finite smooth 2-chain such that dC n Z = 0, 



Proof. We can assume that (7 is a single smooth triangle. Then we can express if on C 
in terms of a complex-valued function 2: and a pointwise-normalized local section (f, i.e. 
if = z(f. Since Im(^) = darg{z), we obtain 

1 diZ 

= -j—pj(dz(/> -I- zV4>,iz(f) = (— 4>,i4>) + {'V(f,i(f) = darg(z) -|- {'V(f,i(f). 

\zf z 

Moreover, away from zeros, we have 

d{V(f, z(f) = {R^(f, i<f) + (V(/) A iV(f) = {R^(f, 1(f) = -H. 

Hence we obtain 

[ = [ daig{z)-\- [ {V<f, i<f) = 27r ^ indexp('0) — [ Q. 

JdC JdC JdC Z 

This proves the claim. □ 

In the case that L is a hermitian line bundle with connection over a closed oriented surface 
M, Theorem 8 tells us that /j^ H G 27rZ. This yields a well-known topological invariant - 
the degree ofL: 

From Theorem 8 we immediately obtain the famous Poincare-Hopf index theorem. 

Theorem 9. Let L ^ M be a smooth hermitian line bundle over a closed oriented surface. 
Then, if if & r(L) is a section with isolated zeros, 

deg(L) = ind^. 
pSM 
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Now, let US consider the discrete case. In general, a section of a discrete vector bundle 
E —7- X with vertex set V is a map i/;: V —)• E such that the following diagram commutes 


E 



TT 


V- 


id 


■V 


i.e. n o Ip = id. As in the smooth case, the space of sections of E is denoted by r(E). 

Now, let L —7- X be a discrete hermitian line bundle with curvature and let ^p G r(L) be 
a nowhere-vanishing section such that 

(8.1) Vijii’i) + -^j 

for each edge ij of X. Here r] denotes the connection of L as usual. The rotation form 
of is then defined as follows: 

Ip 




Remark 10: Equation (8.1) can he interpreted as the condition that no zero lies in the 
1-skeleton ofX (compare Section 11). Actually, given a consistent choice of the argument 
on each oriented edge, we could drop this condition. Figuratively speaking, if a section has 
a zero in the 1-skeleton, then we decide whether we push it to the left or the right face of 
the edge. 

Now we can define the index form of a discrete section: 

Definition 8 . Let L —X 6e a discrete hermitian line bundle with curvature fi. For 
xp G r(L), we define the index form of xp by 

mdA=i(d{* + fi). 

Theorem 10. The index form of a nowhere-vanishing discrete section is h-valued. 

Proof. Let L be a discrete hermitian line bundle with curvature and let t] be its connection. 
Let xp G r(L) be a nowhere-vanishing section. Now, choose a connection form a;, i.e. 
r] = oj(3, where /3 is a trivial connection on L. Then we can write xp with respect to a 
non-vanishing parallel section (p of /3, i.e. there is a C-valued function z such that xp = z(p. 
-) and thus 

i / 


Then ft = arg^^^^^^ 


exp(27rzde4) = exp(z arg(-^) + ^ arg(^) + ^ arg(-^)) = 


^ki^k 


^jk^j ' 


diO. 


'ijk 


Thus 




= 1. 


This proves the claim. 


□ 


If L is a discrete hermitian line bundle with curvature Ll over a closed oriented surface X, 
then we define the degree of L just as in the smooth case: 
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Here we have identified X by the corresponding closed 2-chain. From Theorem 6 we obtain 
the following corollary. 

Corollary 2. The degree of a diserete hermitian line bundle with eurvature is an integer: 

deg(L) G Z. 

The discrete Poincare-Hopf index theorem follows easily from the definitions: 

Theorem 11. Let L —)■ X 6e a discrete hermitian line bundle with curvature Ll over an 
oriented simplicial surface. If tp ^ r(L) is a non-vanishing discrete section, then 

deg(L) = 

ijk^X 

Proof. Since the integral of an exact form over a closed oriented surface vanishes, 

27rdeg(L) = f Q= f d^^ + Ll = 2tt ^ ind^;., 

ijk&X 

as was claimed. □ 


9. Piecewise-Smooth Vector Bundles over Simplicial Complexes 

It is well-known that each abstract simplicial complex X has a geometric realization which 
is unique up to simplicial isomorphism. In particular, each abstract simplex is realized 
as an affine simplex. Moreover, each face a' of a simplex a £ X comes with an affine 
embedding 

ta'a: cr' ^ a. 

In the following, we will not distinguish between the abstract simplicial complex and its 
geometric realization. 

Definition 9. A piecewise-smooth vector bundle E over a simplicial complex X is a topo¬ 
logical vector bundle tt : E —)■ X such that 

a) for each a £X the restriction Eo- := E|^ is a smooth vector bundle over a, 

b) for each face a' of a £ X, the inclusion E^' E^ is a smooth embedding. 

As a simplicial complex, X has no tangent bundle. Nonetheless, differential forms survive 
as collections of smooth differential forms defined on the simplices which are compatible 
in the sense that they agree on common faces: 

Definition 10. Let P he a piecewise-smooth vector bundle over X. An E-valued differen¬ 
tial k-form is a collection oj = {cuo- G Of {a, Eo-)}o-ex such that for each face a' of a simplex 
a £ X the following relation holds: 

where : a' “—)■ a denotes the inclusion. The space of E-valued differential k-forms is 
denoted by fIp^(X, E). 

Remark 11: Note that a 0-form defines a continuous map on the simplicial complex. 
Hence the definition includes functions and, more generally, sections: A piecewise-smooth 
section of E is a continuous section if-.X^E such that for each simplex a £ X the 
restriction V’o- := i/’U • Eo- is smooth, i.e. 

rp,(E) := {f;:X^E\f;^£ r(E^) for all a £ X}. 
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Since the pullback commutes with the wedge-product A and the exterior derivative d of 
real-valued forms, we can define the wedge product and the exterior derivative of piecewise- 
smooth differential forms by applying it componentwise. 

Definition 11. For oj = e V = {^o-lo-eX £ 

U! Ar] := {uJa Ar]a}a£X, du! := {duJa}a£X- 

All the standard properties of A and d also hold in the piecewise-smooth case. 

Definition 12. A connection on a piecewise-smooth vector bundle E over X is a linear 
map V: rps(E) —)■ np^(X, E) such that 

V(/iA) = + / ViA, for all f G n%{X,R), G Tps{E). 

Once we have a connection on a smooth vector bundle we obtain a corresponding exterior 
derivative d^ on E-valued forms. 

Theorem 12. Let E be a piecewise-smooth veetor bundle over X. Then there is a unique 
linear map d^ : Op^(X, E) —)■ np+^(X, E) such that d^fj = Vif for all il) G rps(E) and 

d^ {uj Ar]) = doj Arf -[- { — A d^rj 

for all u G Op^(X,M) and rj G Op^(X, E). 

The curvature tensor survives as a piecewise-smooth End(E)-valued 2-form: 

Definition 13. Let E ^ X be a piecewise-smooth vector bundle. The endomorphism¬ 
valued curvature 2-form of a connection V on E is defined as follows: 

d^ od^ e nl,{X,End{E)). 


10. The Associated Piecewise-Smooth Hermitian Line Bundle 


Let L —7- X be a piecewise-smooth hermitian line bundle with connection V over a sim- 
plicial complex. Just as in the smooth case the endomorphism-valued curvature 2-form 
takes values in the skew-adjoint endomorphisms and hence is given by a piecewise-smooth 
real-valued 2-form Ll: 

d^ od^ = -in. 

Since each simplex of X has an affine structure, we can speak of piecewise-constant forms. 

The goal of this section will be to construct for each discrete hermitian line bundle with 
curvature a piecewise-smooth hermitian line bundle with piecewise-constant curvature 
which in a certain sense naturally contains the discrete bundle. We first prove two lem¬ 
mata. 


Lemma 4. To eaeh closed discrete real-valued k-form oj there corresponds a unique piecewise- 
constant piecewise-smooth k-form oj such that 



for all c G Cfc(X, Z). 


The form oj will be ealled the pieeewise-smooth form associated to oj. 


Proof. It is enough to consider just a single n-simplex a. We denote the space of piecewise- 
constant piecewise-smooth /c-forms on cr by and the space of discrete /c-forms on a by 
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Consider the linear map F: —)■ that assigns to cD G the discrete /c-form 

given by 



Clearly, F is injective. Moreover, since each piecewise-constant piecewise-smooth form is 
closed, we have im F C kerd^, where dk denotes the discrete exterior derivative. Hence it 
is enough to show that the space of closed discrete /c-forms on a is of dimension (^). This 
we do by induction. Clearly, dimkerdo = 1 = (o)- Now suppose that dim kerdfc_i = 

By Lemma 2, we have kerd^ = imdfc_i. Hence 

dimkerdfc = dimimdfc_i = dim 12^ - dimkerdfc_i = ("'^^) - = Q). 

Therefore, for each closed discrete fc-form we obtain a unique piecewise-constant piecewise- 
smooth /c-form which has the desired integrals on the /c-simplices. □ 

It is a classical result that on star-shaped domains U C each closed form is exact: If 
n G M) is closed, then there exists a form a; G M) such that = doo. More¬ 
over, the potential can be constructed explicitly by the map K: 17^(17, M) —)■ M) 

given by 

k _ 

K{^) = ^ t''~^^ii-i^,{tx)dt^ Xi^dxi^ A ... Adxi^ A ... Adxi^, 

where 17 = X]*i<---<4 ^h-ik A ... A dxi^. One directly can check that 

K{d9.) + dK{yL) = n. 

Hence, if dl7 = 0, we get 17 = dK{Q). The same construction works for piecewise-smooth 
forms defined on the star of a simplex. This yields the following piecewise-smooth version 
of the Poincare-Lemma. 

Lemma 5. On the star of a simplex each closed piecewise-smooth form is exact. 

Now we are ready to prove the main result of this section. 

Theorem 13. Let L —)■ X 6e a discrete hermitian line bundle with curvature 17 over a 
simplicial complex and let 17 he the piecewise-smooth 2-form associated to 17. Then there 
is a piecewise-smooth hermitian line bundle L —)■ X with connection V of curvature 17 such 
that Li = Li for each vertex i and the parallel transports coincide along each edge path. 
The bundle L is unique up to isomorphism. 

Proof. First we construct the piecewise-smooth hermitian line bundle. Let L —)■ X be a 
discrete hermitian line bundle with curvature 17 and let ij denote its connection. Let V be 
the vertex set of X and let Si denote the open vertex star of the vertex i. Further, since 17 
is closed, by Lemma 4, there is a piecewise-constant piecewise-smooth form 17 associated 
to 17. Now, consider the set 

L := U Si X Li. 
i&V 

Note, that Si n 5^ / 0 if and only if ij is an edge of X or i = j. Thus, if we set iqu := id|L., 
we can define an equivalence relation on L as follows: 

ii,p,u) ~ ij,q,v) p = q and v = exp(-z / ^ n)r]ij{u), 

^3 
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where Afj denotes the oriented triangle spanned by the point i, j and p. Note here that A^j 
is completely contained in some simplex of X. Let us check shortly that this really defines 
an equivalence relation. Here the only non-trivial property is transitivity. Therefore, 
let {i,p,u) ~ {j,q,v) and {j,q,v) ~ {k,r,w). Thus we have p = q = r and p lies in a 
simplex which contains the oriented triangle ijk. Clearly, the 2-chain AC + -|- A^- is 

homologous to ijk and since piecewise-constant forms are closed we get 


H = - 


/A^’.+AL 

ij jk 


I kl f — j H + Q-ijk- 
JaI. Jijk ./al 

k.i. '' i.k. 


Hence we obtain 


w 


= exp(-z J ^ Q.)rijk(exp[-i j ^ ^)rjij{u) 

^3 k ij 

= exp(-z^^ ^)r]jkor]ij{u) 

jfe 

= exp(-z J ^ ^ - iVLijk)rijk o riij{u) 

‘~^ik 

= exp(-z / ^ n)pik{u), 


and thus {i,p,u) ~ {k,r,w). Hence ~ defines an equivalence relation. One can check 
now that the quotient L := \jj^ is a piecewise-smooth line bundle over X. The local 
trivializations are then basically given by the inclusions 5* x Lj L sending a point to 
the corresponding equivalence class. Moreover, all transition maps are unitary so that the 
hermitian metric of L extends to L and turns L into a hermitian line bundle. Clearly, 



= L. 


Next, we need to construct the connection. Therefore we will use an explicit system of local 
sections: Choose for each vertex i G V a unit vector Xi G Lj and define (j)i{p) := [i,p,Xi]. 
This yields for each vertex i a piecewise-smooth section cji define on the star Si. For each 
non-empty intersection Si n 5^ / 0 we then obtain a function gij ; Si H Sj —)■ S. By the 
above construction, we find that, if rjij^Xi) = VijXj, 


( 10 . 1 ) 


9ij{p) 


Vij exp 



Since kl is closed. Lemma 5 tells us that ^\si is exact. Hence there is a piecewise-smooth 
1-form LJi defined on Si such that dui = fils'.. In general, the form uji is only unique up 
to addition of an exact 1-form, but among those there is a unique form uJi which is zero 
along the radial directions originating from i. To see this, just choose some potential Ai 
of fl|sj and define a function /: 5* —)• M as follows; 

For p G Si, let f{p) := Lp oji, where y/’ denote the linear path from the vertex i to the 

'i 

point p. Then oji := Wj — df is a piecewise-smooth potential of and vanishes on radial 
directions. For the uniqueness, let oJi be another such potential. Then, the difference 
oji — uJi is closed and hence exact on Si, i.e. there is /: 5* —)• M such that df = uJi — uJi. 
Since df vanishes on radial directions / is constant on radial lines starting at i and hence 
constant on Si. Thus Ui = uii. 


Suppose that for each edge ij the forms uJi and ujj are compatible, i.e., wherever both are 
defined. 


lUJj = ICOi + d log Qij . 
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Then we can define a connection V as follows: Let i/: G and let X G TpCJ for some 

simplex a of X, then there is some Si 3 p. On Si we can express ip with respect to (pi, i.e. 
Ip = z (pi for some piecewise-smooth function z: S'j —)• C. Then we define 

(10.2) VxV’ := {dz{X) - iu}i{X)z)(pi. 

In general there are several stars that contain the point p. From compatibility easily 
follows that the definition does not depend on the choice of the vertex. Hence we have 
constructed a piecewise smooth connection V. One easily checks that V is unitary and 
since doJi = Ills, we get d^ o = —zQ, as desired. 

So it is left to check the compatibility of the forms Uij constructed above. Let ij be 
some edge and let po be a point in its interior. Since Ui — ojj is closed, we can define 
ip: n S'j —)■ M by (/9(p) := oji — uj, where 7 p is some path in Si O Sj from the point po 

to the point p. Then, for p G SiCl Sj, 


Q = I CJ, = / LOj = 

Jij+lj-lp 


Id A, 


- = / UJi- UJj = ip{p), 


where as above denotes the linear path from i to p and, similarly, 7? denotes the linear 
path from j to the point p. From this we obtain 


oJi — (jjj = dip = d Q 

J 

and in particular zujj = zuzi + dloggij. This shows the existence. 

Now suppose there are two such piecewise-smooth bundles L and L with connection V 
and V, respectively. We want to construct an isomorphism between L and L. Therefore 
we again use local systems. Explicitly, we choose a discrete direction field X G L. This 
yields for each vertex i a vector W G Lj = Lj which extends by parallel transport along 
rays starting at f to a local sections (pi of L and, similarly, to a local section (pi of L defined 
on Si- 

Now we define F: L —)■ L to be unique map which is linear on the fibers and satisfies 
F{<pi) = (pi on Si- To see that F is well-defined, we need to check that it is compatible 
with the transition maps. But by construction both systems have equal transition maps, 
namely the the functions gij from Equation (10.1) with Vij given by rnj^Xi) = VijXj. Now, 
if Zi (pi = Zj (pj , then Zi = Zj gij and hence 


F{zipi) = Zi (pi = Zigijpj = Zj (pj = F{zj pj). 


Using Equation (10.2) one similarly shows that F o V = V o F. Thus L = L. 


□ 


11. Finite Elements for Hermitian Line Bundles With Curvature 

In this section we want to present a specific finite element space on the associated piecewise- 
smooth hermitian line bundle of a discrete hermitian line with curvature. They are 
constructed from the local systems that played such a prominent role in the proof of 
Theorem 13 and the usual piecewise-linear hat function. 

Let L be the associated piecewise-smooth bundle of a discrete hermitian line bundle L —)■ X 
and let Xj: X —7- M denote the barycentric coordinate of the vertex i G V, i.e. the unique 
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piecewise-linear function such that Xi{j) = 6ij, where 6 is the Kronecker delta. Clearly, 


r(L) = 0L,. 

ieV 


To each X £ hi we now construct a piecewise-smooth section if: as follows: First, we 
extend X to the vertex star Si of the vertex i using the parallel transport along rays 
starting at i. To get a global section V’ £ rp 5 (L) we use Xi to scale 4> down to zero on dSi 
and extend it by zero to X, i.e. 


i’p 



forp £ Si, 
else. 


The above construction yields a linear map l: r(L) (L). Clearly, t is injective - a 

left-inverse is just given by the restriction map 


rps(L) 3 'tjj G r(L). 

Definition 14. The space of piecewise-linear sections is given by rp/(L) := im/.. 


Thus we identified each section of a discrete hermitian line bundle with curvature with a 
piecewise-linear section of the associated piecewise-smooth bundle. This allows to define 
a discrete hermitian inner product and a discrete Dirichlet energy on r(L), which is a 
generalization of the well-known cotangent Laplace operator for discrete functions on tri¬ 
angulated surfaces. Before we come to the Dirichlet energy, we define Euclidean simplicial 
complexes. 

Similarly to piecewise-smooth forms we can define piecewise-smooth (contravariant) k- 
tensors as collections of compatible A:-tensors; A piecewise-smooth k-tensor is a collection 
T = {To-jo-ex of smooth contravariant fc-tensors T„ on a such that 

= T,,, 

whenever cr' is a face of a. A Riemannian simplicial complex is then a simplicial complex 
X equipped with a piecewise-smooth Riemannian metric, i.e. a piecewise-smooth positive- 
definite symmetric 2-tensor g on X. 

The following lemma tells us that the space of constant piecewise-smooth symmetric 
tensors is isomorphic to functions on 1 -simplices. 

Lemma 6. Let X be a simplicial complex and let £ denote the set of its 1-simplices. 
For each function /:£—)• M there exists a unique constant piecewise-smooth symmetric 
2-tensor S such that for each 1-simplex e = {i,j} 

Se{j -i,j -i) = /(e). 


Proof. It is enough to consider a single affine n-simplex a = {io,..., in} with vector space 
V. Consider the map F that sends a symmetric 2-tensor S' on 1/ to the function given by 

F{S){e) := S{ik - ij,ik - ij), e = {ij,ik} C a. 

Clearly, F is linear. Moreover, if Q denotes the quadratic form corresponding to S, i.e. 
Q{X) := S(A, A), then 

Six, Y) = i(Q(A) + QiY) - QiX - Y)). 
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Hence, from F{S) = 0 follows S = 0. Thus F is injective. Clearly, the space of symmetric 
bilinear forms is of dimension n{n + l)/2, which equals the number of 1-simplices. Thus 
F is an isomorphism. This proves the claim. □ 

It is also easy to write down the corresponding symmetric tensor in coordinates; Let 
a = {io,..., in} be a simplex. The vectors ej := ij — zq, j = 1,... ,n, then yield a basis 
of the corresponding vector space. Let / be a function dehned on the unoriented edges of 
a and let Xi- denote the barycentric coordinates of its vertices ij, then the corresponding 
symmetric bilinear form S[ is given by 

(11.1) Sl= fioij dXi^ ^ dXi. + fioik- fijik)dXi^^dXi^. 

Thus starting with a positive function /, by Sylvester’s criterion, it has to satisfy on each 
n-simplex n — 1 inequalities to determine a positive-definite form. If the corresponding 
piecewise-smooth form is positive-definite, we call / a discrete metric. 

Definition 15. A Euclidean simplicial complex is a simplicial complex X equipped with 
a discrete metric, i.e. a map i that assigns to each 1-simplex e a length £e > 0 such that 
for each simplex a the symmetric tensor is positive-definite. 

Now, let X be a Euclidean simplicial manifold of dimension n and denote by X^ the set of 
its top-dimensional simplices. Since each simplex of X is equipped with a scalar product 
it comes with a corresponding density and hence we know how to integrate functions over 
the simplices of X. Now, we define the integral over X as follows: 



Moreover, given a piecewise-smooth hermitian line bundle L —)■ X with curvature, then 
there is a canonical hermitian product ((.,.)) on rps(L): If i/;, 0 G rps(L), then 

= [ {'4’A)- 

Jx 

In particular, if L is the associated piecewise-smooth bundle of a discrete hermitian line 
bundle L with curvature 11, then we can use i to pull ((.,.)) back to r(L). Since i is injective 
this yields a hermitian product on r(L). 

Now we want to compute this metric explicitly in terms of given discrete data. 

Definition 16. A piecewise-linear section if G rpj(L) is called concentrated at a vertex i, 
if it is of the form = t(V’i) for some vector ifi Ghi. 

It is basically enough to compute the product of two such concentrated sections. Therefore, 
let fii G Lj and G Lj and let and denote the corresponding piecewise-linear 
concentrated sections. 

Now consider their product Clearly, this product has support SiCiSj. For 

simplicity, we extend the discrete connection r] to arbitrary pairs ij in such way that r]ii = 
id and rjij: Lj —)■ Lj is zero whenever {i,j} 0 X. With this convention. Equation (10.1) 
yields 

= {i>j,r}ij{'fii)) XiXj exp(-^ / ^ h), 

J . 


( 11 . 2 ) 
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where Q denotes the constant piecewise-smooth curvature form associated to fi. 

Now, let us express the integral over A^- on a given n-simplex. Therefore consider an 
n-simplex a = {io, • • • ,in}- The hat functions Xjj,... ,Xi^ yield affine coordinates on a 
and we can express any 2-form with respect to the basis forms dxi- A dxi^. One can show 
that 


UJ.J.CAIU 

/ dxi. A dxi^ = 

J a' 


for a' = ±ijikie, 

Thus we obtain 

1 0 

else. 


n= Y, 

2 ^ dXij^ . 


l<j<k<n 

Now we want to compute the integral over the triangle A^^^^ C a. By Stokes theorem, 

ril rp rio 

dxijAdxi^= / Xi.dxi^+ / Xi.dxi^+ / Xi.dxi^, 

JiQ Jii Jp 

where the integrals are computed along straight lines. A small computation shows 

J ^ dxi^ A ^ (5y X 4 (p) - 6-ik Xi. (p)), 

ion 

Thus, for j < k, we get f^p dxi^ A dxi^ = Idy Xi^{p) and hence 

»on 






dxi. A dxi^ 


^iohijXij {p), 

j 


where we have used the convention that 0 vanishes on all triples not representing an 
oriented 2-simplex of X. With this convention Equation (11.2) becomes 

(11.3) = {'tpj,riij{ipi)) XiXj exp[-i'^nijkXk). 

k 


In particular, using Equation (11.3), we can compute the norm of a piecewise-linear section 
on a given triangle ijk. Therefore we distinguish one of its vertices, say i, and write 
with respect to a section which is radially parallel with respect to i. Now, one checks that 

1^1 = |q + - a) -k - a)], 

where Cj, cj, Cfc G C are constants depending on the explicit form of i/^. An example of the 
norm of a piecewise-linear section is shown in Figure 7. 

As the next proposition shows, the identification of discrete and piecewise-linear sections 
perfectly fits together with the definitions in Section 8 . 

Proposition 3. Let ip G r(L) be a discrete section and let ip G rp;(L) be the corresponding 
piecewise-linear section, i.e. ip = l{iP). Then, if ip has no zeros on edges, the discrete 
rotation form and the piecewise-smooth rotation form are related as follows: For 
each oriented edge ij, 

4 = / A. 

Jij 


Proof. The claim follows easily by expressing ip with respect to some non-vanishing parallel 
section along the edge ij. □ 
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Figure 7. The graph of the norm of a piecewise-linear section of a bundle 
over a torus consisting of two triangles. Its two smooth parts fit contin¬ 
uously together along the diagonal. In this example the curvature of the 
bundle over each triangles is equal to dvr. Note that the section has 4 zeros 
- just as predicted. 


In particular, by Theorem 8, the index form of a non-vanishing section of a discrete her- 
mitian line bundle with curvature counts the number of (signed) zeros of the corresponding 
piecewise-linear section of the associated piecewise-smooth bundle. 

Let us continue with the computation of the metric on r(L). To write down the formula 
we give the following definition. 

Definition 17. Let X be an n-dimensional simplicial manifold and let Q G M). To 

an n-simplex a and vertices i,j,k,l ofX we assign the value 




vol I 


cr 


XkXi exp -z 






i), 


where have chosen for integration an arbitrary discrete metric on X. 

Remark 12: Note that the functions indeed well-defined. On a simplex, any 

two such measures induced by a discrete metric differ just by a constant. 

With Definition 17 and Equation (11.3) we obtain the following form of the metric: 

Theorem 14 (Product of Discrete Sections). Let h be a discrete hermitian line bundle 
with curvature Ll over an n-dimensional Euclidean simplicial manifold X, then the product 
on r(L) induced by the associated piecewise-smooth hermitian line bundle is given as 
follows: Given two discrete sections i/jj, f fi, 

{{'•PA)) = -where An = ©^AjCbl) vol (a). 

{i,j}CcreXn 

Note that Q^^ j{k,l), and hence AA computed explicitly using Fubini’s theorem 

and the following small lemma, which can be shown by induction. 
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Lemma 7 . Let c G C*, n G N and [a, 6 ] C M 6 e an interval. Then 

x-exp(cx)dx = exp(cx) 

Next, we would like to compute the Dirichlet energy of a section xjj G rpj(L), i.e. 

EdW = [ 

Jx 

Note, that the Dirichlet energy comes with a corresponding positive-semidehnite hermitian 
form ((., .))£) - called the Dirichlet product. Clearly, like the metric, the Dirichlet product 
is completely determined by the values it takes on concentrated sections. 

In general, if 7 ^ 0 is piecewise-linear section concentrated at i, it is given on the vertex 
star Si as a product = Xi <p, where Xi denotes the barycentric coordinate of the vertex 
i and i;A is a local section radially parallel with respect to i. Clearly, 

= dxi ^ + iXiOJi^, 

where w* denotes the rotation form of (j), i.e. Vcp = iuji4>. Note here that Ui does not 
depend on the actual value of ■0 at i, but is the same for all non-vanishing piecewise-linear 
sections concentrated at i. 

To compute the rotation form w* at a given point po G Si, we use a local section ( which 
is radially parallel with respect to po such that Cpo = 0po • Then we can express 0 in terms 
of (, i.e. 

0 = zC, 

for some piecewise-smooth C*-valued function z dehned locally at po. Clearly, \z\ is 
constant, and hence 

ho = = Cpo +^(Po) VCIpQ = dlogz\p^ 4 >po = idargz\p^ 0 po. 

The clue is that we can now use the relation of parallel transport and curvature to obtain 
an explicit formula for z. If p is sufficiently close to po, then the three points p, i and 
Po determine an oriented triangle which is contained in a simplex of X. Its boundary 
curve 7p consists of three line segments 71,72,73 connecting p to i, i to po and po back to 
p. Hence on each of these segments either 0 of 0 is parallel and 



Thus we obtain that z{p) = exp(—z f^p D) and hence 

LUi = —/ D) 
po po 

Now, if is contained in a simplex a = {zq, . •., z„}, one verifies that 

dxi. A = ^{xi^{po)xhiP) - Xi^:{po)xij{p)). 
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Thus, 


d(^ / ^ / dxij A dxi^,^ 

PO l<j<k<n 

= H ^ioijik(xi^dxi^-Xi^dxi.) 
l<j<k<n 


PO 


PO 


— '^2 (y^ ^iQiiik^ik)dXii 

l<j<n k^j PO 

and, using the convention on 0, from above, we find the following simple formula: 
(11.4) CUj — ^ ) ( y~) Tfc) dXj 

j k 

where we sum over the whole vertex set of X. 


Now, given this local form expressions, we can finally return to the computation of the 
products which we are actually interested in. Therefore we consider two piecewise-linear 
sections concentrated at the vertices i and j: 

y* := := 

for some ipi G Lj and 'ijjj G Lj. On their common support Si O Sj both section can be 
expressed, just as above, as products of a real-valued piecewise-linear hat functions Xi and 
Xj and radially parallel local sections cj)^ and (fS: 

■0* = Xj (^*, = Xj ^. 


Clearly, 


= / {dxjcjS -|- iXjUOj (j)^ , dxicj)^ + iXiiOi (j)^) 
JSi n Sj 


= / {dxj + iXjijjj, dxi + iXiOJi) {(jS, (//). 
JSi n Sj 

With Equation (11.3) we see that 

= (V’i,%(V'i)) exp(-z E ^ijmXr, 


Moreover, by Equation (11.4), 


{dxj + iXjOJj, dxi + iXiOJi) 


+ i\2^{VLik'i'XiXii{dxj,dxk') - VLjk'i'XjXv{dxk'-,dxi)) 
k',V 


The constants {dxk',dxii) are basically provided by the following lemma. 

Lemma 8. Let a = {vq, ... ,Vn} be a Euclidean simplex of dimension n > 0 and let Xj 
denote its barycentric coordinate functions. Then 

gradxj = -^Ni, 
hi 

where hi denotes the distance between Vi and ai = a \ {xj} and Ni denotes the outward¬ 
pointing unit normal of ai. 
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Proof. This immediately follows from two basic facts; First, dxi{vj — vq) = Sij for i,j > 0. 
Second, hi = {vq — Vi, Ni). □ 


Lemma 8 yields almost immediately a higher dimensional analogue of the well-known 
cotangent formula for surfaces. 

Theorem 15 (Cotangent Formula). Let a be a simplex of a Euclidean simplicial complex 
X and let dimer >1. If i j, 


c 


ij 

<7 


{dxi, dxj) 


n{n—l) 


cot a*^vol ((T\{i,j}), 


if {hi] C cr, 
else. 


Here denotes the angle between the faces a \ {i} and a \ {j}. Moreover, 




PI := / \dxif = , 

Ja (0 else, 

where hi denotes the distance between the vertex i and the face a \ {i}. 


Proof. Clearly, if {i,j} (f o, then f^{dxi,dxj) = 0. Now, let {i,j} C a, i j. With the 
notation of Lemma 8, we have 



(grad Xi, grad Xj) vol a 


{Ni,N,) 

hihj 


vol a. 


Furthermore, cosa^ = —{Ni,Nj) and n!volcj = (n — 2)!sina^ vol (cr \ {i, j}). This 
yields the first part of the theorem. Similarly, nvola = hi vol (cr\ {f}). Setting then i = j 
yields the second part. □ 


Definition 18. Let X be an n-dimensional simplicial manifold and let G n^(X, M). Let 
a be an n-simplex and i,j,k,l be vertices ofX. Then, let 


. . ■= 




w \ /* GXpf ^ ^ ^ijm^m) ? 

vol (a) J a ^ 

vol (a) I exp(-z^ny>„Xm), 


where we choose for the integration an arbitrary discrete metric on X. 

Remark 13: Just like the functions ike values and the functions and 

are well-defined (compare Remark 12). 

Now, with these definitions, we can summarize the above discussion by the following 
theorem. 

Theorem 16 (Discrete Dirichlet Energy). Let L be a discrete hermitian line bundle with 
curvature D over an n-dimensional Euclidean simplicial manifold X, then the Dirichlet 
product on F(L) induced by the associated piecewise-smooth hermitian line bundle is given 
as follows: If 4> = J2i '/’* V’i iwo discrete sections, 

((</>, V’))d = ^ ^VijAi)), w'H = ’ 

hi {i,i}C(TeX„ 
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where 

(11.5) 


■ = 


+ 4"' ©)?,,,(i, I') - n.k'v cf I')) 

k',i' 


12. Discrete Energies on Surfaces - An Example 


While the computation of the Dirichlet product {{■, ■))d and the metric ((.,.)) of discrete 
sections is quite complicated and tedious for higher dimensional simplicial manifolds, it is 
manageable for the 2-dimensional case. We are going to compute it explicitly. 

Throughout this section let L denote a discrete hermitian line bundle with curvature D 
over a Euclidean simplicial surface X and let a = {i,j, k} be one of its triangles. 

The metric ((.,.)) is easily obtained. We basically just need to compute the values 0^ j j(i, i) 
and which can be done over the standard triangle. We get 


( 12 . 1 ) 


0 


n 


.(b^) = g, 


0 ?«(bj) = 2 


exp( 1 -\- -j- 


iQ 2 _ Dq3 




Now, we compute the Dirichlet product ((., .))d on X. For n = 2, the expressions ^ 
and simplify drastically. First, we look at the diagonal terms. We have 


E 

k',k\i',i’ 


cE ^ik'l'^ikH" 






= k) - 2ci^E^,,(j, k) + cf 


and with 

we get the following formula: 


W^- ■ = 

' G.l.l 


+ 


rJJ 


— -L 

^(T ' ^tj f)2 

-Djo “b'fc- 


Now we would like to obtain a similar formula for the off-diagonal terms. Since dxi+dxj = 
—dxk, we have Hence, 


E ^ ^ik'l'^jk'l' I') 

k',k’,1',1’ 





{k,k) + j) + r. 






This time the expressions become more complicated. We get 

- SDf-fc \i^%k + (“20 - ^i^ijk + exp(-i%fc)), 

ijk ^ ' 

2?,i,i(uj) +S^,ij(Tfc) = ^(-6-PTflyfc -t Xyj4 _ Tzflf.fc -t (6-t exp(-zflyfc)). 

ijk ^ ' 
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Thus, 


-|-([6c“ - 20c«] + [I2c« - + [3c« - - f - f fij, 

^^ijk 

+ C+ ([20c^^ - 6c“] + [8c*^ - exp(-z%fc). 

Now, let us look at the second sum in Equation (11.5). We have 


k',l' 




The formula for Ahj) is already given in Equation (12.1). Eurther, we have 


o-,*J 

2 


*:) = 04 “ (3 - 2!S2yid - + (-3 + ifiyt) exp{-sSiyt)) = i). 




Thus we get 

©?,,,(j,r)) = 


fc',Z' 

2 


fit 




([3(c" + ciJ) - c“].f!yt + [2(c“ + <i>) - c“]SjJ, + 1 [cf - c" - 


+ ^-S2h + ([c“ - 3(c; + ci’)]tSJyt+[c; +<!i']Sjj 4 .)exp(-!S!yj:)). 


Hence, with 


^cr,i,j ~ q 4 {^ijk '^^ijk flijfc ®^p( 


Equation (11.5) becomes 


^ ( [6<"“ - 20ct] + [l2c(f' + 3(c" + c") - 5cf ] + [4c(?' + 2(c“ + c>J - cf)] Hi 

ijk 


2 

ijfc 


+ 1 [3(cf - c" - ci^ ) - 84^] 


Ufe 


+ ( - 6c“] + - 3{c“ +c’J)- c^''] + [c“ - 2c^ + cj-^] 12^^,) exp(-*ny fc) ). 

Since n = 2, the weights are just given as follows: 

ij _ cot kk _ ^ij 

2 ’ " 2/ifc’ 

where iij denotes the edge length. We would like to express them explicitly in terms of 
the Euclidean metric of a. In fact, we can distinguish the vertex k as origin and use 
the hat functions Xi and Xj as coordinates on a. With respect to these coordinates, the 
metric is given by a matrix: 

9 = 


9ii 912 
921 922 
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In terms of g the cotangent weights are given as follows: 


512 


= 


2^Jdeig'’ 

9 ii - 2512 + 522 


(jk = - 


2-y/det g ’ 
and we have rederived the formulas in 


= 


511 - 512 

2-v/det g 
522 


^ki ^ _ 522 - 512 


2\/det g ’ 


ci’ = 


2v^det5 ’ 
511 


2-y/det (/’ 


■ =- 
vol ( 


( [3511 + 4512 + 3^22] — [511 + 512 + 522] l^ijk + 

+ - ([3511 + 4512 + 3^22] 

+ [2511 + 3c/i2 + 2 g 22 ] t^ijk - 5 [511 + 2512 + 522 ]f^|fc) exp(-^^l^Jfc))• 
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